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Abstract. To model the decay of a quasibound state we use the modified two-potential 
approach introduced by Gurvitz and Kalbermann [1,2]. This method has proved itself 
useful in the past for calculating the decay width and the energy shift of an isolated 
quasistationary state [5]. We follow the same approach in order to propagate the 
wave-packet in time with the ultimate goal of extracting the momentum-distribution of 
emitted particles. The advantage of the method is that it provides the time-dependent 
O | wave function in a simple semi-analytic form. We intend to apply this method to the 

modeling of metastable states for which no direct integration of the time-dependent 
Schrodinger equation is available today. 
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The Two Potential Approximation (TPA) introduced in Refs. [1,2] turned out to 
be an extremely successful tool for the description of a metastable state. Simple 
expressions based on the TPA made it possible to obtain a very precise estimate of 
the life-time of a very narrow resonance without the need of introducing an explicit 
time dependence. In this work, we use the TPA in order to derive the equations 
describing the time evolution of the wave function of a particle tunneling through 
a spherically-symmetric barrier. 

Let us consider a particle moving in a central potential V(r) with a barrier. 
Asymptotically, i.e., at large values of r, we assume that V{r)— >0. In the TPA, 
V(r) can be decomposed as 

V(r) = U(r) + W(r), (1) 
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where 



u(r) - { V (r) Xr<R 

U[r) ~ \ V(R) \lr>R [Z) 

is an auxiliary potential that produces a bound state at energy Eq close to the 
energy of the metastable state, and W(r) is a "closing" potential which is treated 
perturbatively. The separation radius R should be chosen far from the classical 
turning points [3]. 

At t—0 the initial state is taken to be the bound eigenstate ®o(r) of the auxiliary 
Hamiltonian 



H = T + U(r) (3) 

(we take h—1). In the following we assume that <&o(r) is well isolated, i.e., it is 
well separated from the remaining bound states of U(r) having the same quantum 
numbers. In such a case, at t>0, the wave packet represented by the wave function 
\l/(f, t) can be expanded in the basis {<3> (r), $ k (r)}: 

*(?,t) = b (t)M^- iEot + J b k (t)^ k (f)e- iE ^ , (4) 

with the initial conditions bo(t = 0)=1 and b k {t = 0)=0. In Eq. (4) the wave 
functions $fe(r) represent the continuum and E k = V(R) + k 2 /2m. We shall refer 
to the first and second terms above, as ^/(r, t) and ^//(r, t), respectively. 

To evaluate the two components, Wj(f, t) and ^//(r, t), the Laplace transform 
method can be applied. In terms of the Laplace-transformed expansion coeffi- 
cients bit), 

roo 

6(e) = / b(t)e l£t dt, (5) 
Jo 

the Laplace transform of the wave packet \I/(r, t) can be written as 

y I (r,e + E )=b (e)%(r), (6) 
* n (r,e + E )=J^- 3 l k (e k ) $ fc (f) , (7) 

where b k (t) = e- iV{R H k (t) and e k = e + E + V(R) - E k . 

Assuming a spherically symmetric potential V(r), the coefficient b (e) has been 
calculated as [2] 



bo(e) = — - — , 
e -e 



(8) 



with 



e = A-i^ (9) 



where a = yj2m(Vo — E ) and A; = y 2m(E + e )- 111 this work, 0o( r ) is the 

radial wave function of ^ an d X«fe( r ) an d XzA^( r ) are > respectively, the regular and 
outgoing waves of the Hamiltonian with the potential W(r) = W(r)+V(R). (Note, 
that our radial continuum functions satisfy the orthogonality and completeness 
relationships 

roc 

/ Xlk(r)xw(r) dr = 5(k - k') , (10) 
Jo 

/ X*i k {r)xik{r')dk = 6(r-r'). (11) 
J o 

Compared with expressions in Refs. [1,2], this results in an additional factor of 
\Jk/2 in the front of every Xik(j') function [4].) 
With the above definitions, the radial part of the first component in Eq. (4) is 

^( r> t) = Mle^f. (12) 
r 

The coefficients b k (e k ) are determined by solving the system of integral equations 

f d 3 k ~ - 

eb (e) = % + W m b^{e) + j W ok b k (s k ) , (13) 

f d 3 k' ~ - 
e k b k (e k ) = W ko b (e) + y ^— ^ W fcfc /6fc/(e fc /) , 

with W^fe/ = {<fr k \W\& k >). The solution of (13) can be formally written as 
h(s k ) = — <$ fc | (l + + W^o^Go + • • • ) W | $ > 6o(e) , 

where 

G °-i(2^F^^- (14) 

The outgoing part of the wave function, ^n(f, e), can now be expressed in terms 
of */(r, 5) as 

* 77 (f, £ + E ) = J d 3 r' G(e + E ; 7, 7) W(r') ^(f, e + E ) , (15) 



where we now introduce the Green's function 

G(E) =G (E) + Gq(E)WG(E) = (I - k)(E - H + KW)-\ (16) 

with A = |$ ) (<&o| being the projection operator on <3> - The Green's function G(E) 
is approximated in the spirit of Ref. [2] by neglecting the contribution from A, and 
then by replacing the potential V(r) by W(r). This gives G(E) w Gy / (E) ) where 

G W {E) = {E - H w y l (17) 

is the Green's function of = T + W. 

By taking the inverse Laplace transform of (7), one obtains for the the radial 
wave function 

| poo rt7+oo _ 

^ n (r, t) = — r'dr' W{r')^{r') / de e~ t£t G w (e; r, r')b (e - E ) . (18) 

The e-integral is evaluated using the residue theorem, and results in the sum of the 
residues corresponding to the two poles of the integrand. 

Contribution due to the pole of b (e — E ) 

Using the standard techniques explained in Ref. [2], we obtain 

Ma(r < R,t) = y| [a^(R) + dg'(R)] Xlk0 (r) e^**** , (19) 



and 



+ vf ^^o(R) + WtfW^ • (20) 

Note that for r > R the contribution from t/jj is exactly canceled by the first term 
in (20). 

Contribution due to the pole of the Green's function. 

The Green's function G^ has a continuum of simple poles along the real E > 
axis. After using the spectral representation of G^r{e;r, r'), one can express 
ipn,b(r,t) as 

1>iijb(r,t) = — / dr'W(r')Mr') / dk — — - 2 X i*(r) Xi k (r') ■ (21) 

T J R JO rC — /Cq 



The evaluation of the integral (21) represents the corner stone of the present ap- 
proach. 

For now we will restrict ourselves to making some remarks regarding the asymp- 
totic behavior of this integral at large values of r. We shall also assume that the 
potential V(r) has finite range (i.e., it vanishes at large values of r). While this as- 
sumption cannot be used for the case of the Coulomb potential, it is still interesting 
to investigate the general structure of the solution for the short-range potential. In 
this limit, the ^-matrix is meromorphic in the complete complex A;-plane, and [4,6] 

[S l (k)]* = S l (-k*), and Si(fc) = [S^-k)}- 1 . (22) 

Expressing X/fc( r ) m the asymptotic form: 

Mr) = \ % {Sl /2 (k)e^ - (-) l S[ 1/2 (k)e- ikr ) , (23) 
one obtains for the /c-integral in (21): 



-iJs-t 

POO p 1 2m 

I(r, r', t) = j o dk Xlk {r) X * lk {r>) (24) 



fc 2 1 

1 f°° ^ e 2m 

oo kP 1 /lq 



2tt J- 



ik(r-r') _ s^k) e tk(r+r,) . (25) 



The integrand in Eq. (25) is a sum involving two complex functions of complex k, 

1 and #%, (26) 



h.1 _ 1.2 ' h.2 _ p ' 

which have common poles at ±k . In addition, Si(k) has an infinite number of 
simple poles. They are located in the lower half of the complex A;-plane, symmet- 
rically with respect to the imaginary axis. Following the notation of van Dijk and 
Nogami [7], we shall denote the poles in the fourth quadrant with k u , v — 1, 2, 3, . . ., 
and the poles in the third quadrant with k u , v = —1,-2,-3,... It follows from 
Eq. (22) that 

Re(fc„) = -Re{k- U ), Im(k u ) = hn(k- v ) . (27) 

In the following, the residue of the Si(k) at the pole k u is denoted by b u . 

Since the complex function (26) has no essential singularity at infinity, we can 
apply the Mittag-Lefner theorem in order to obtain a pole expansion for (26). 
Consequently, Eq. (26) can be replaced by 

1/1 1 \ 

and (28) 



2ko V k — ko k + ko 
Si(fco) 1 Sti-ko) 1 



oo 



2k k- k Q 2k k + k ^fr^ k 2 v -kl k-k v 



+ U 2 _ 1.2 u u ■ ^ 



By substituting Eqs. (28) and (29) in (25), the integral (25) becomes 

1 1 



I(r, r',t) 



-00 4rrk \k — k k + k 



(") 
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-00 



2k \ k — k k + ko 



dk ' 



E 



U2 _ U2 

u=-oo K v H 



The above can be now expressed in terms of the Moshinsky function 



M(k,TZ,r) = — dp 

Z7T 



roc e~ ip2r e~ ipn 



2-k J -00 1 p — k 



where 



y = e 



(30) 



(31) 



where r = t/2m and 1Z — r ± r'. The integral I(r,r',t) can now be calculated in 
the closed form: 



I(r,r',t) 



2ik 



M (*»•''-'''• i) + M (^ r '- r -i) 

)'S,(k )M (k ,r + r', J-) - {-)'S,(ko)M (k , -r - r', ^) 



(32) 



This concludes our derivation. 
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